Numerical solution of a Model for turbulent Diffusion

A model is considered for turbulent diffusion which consists of a Riesz space fractional derivative to describe the turbulent phenomenon and also includes advection and classical diffusion. We present a first order explicit numerical method and a second order implicit numerical method to solve our problem and prove convergence results for both methods, including the derivation of stability constraints needed for the explicit numerical method to converge. In the end, to give some insights into the phenomenon of turbulent diffusion described by the Riesz fractional derivative, we show the behavior of the solution when we consider a Gaussian initial condition.